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Abstract 

The five-band Hubbard model for a d band with one electron per site is a model which has very 
interesting properties when the relevant ions are located at sites with high (e. g. cubic) symmetry. 
In that case, if the crystal field splitting is large one may consider excitations confined to the 
lowest threefold degenerate t2g orbital states. When the electron hopping matrix element (t) is 
much smaller than the on-site Coulomb interaction energy (U), the Hubbard model can be mapped 
onto the well-known effective Hamiltonian (at order /U) derived by Kugel and Khomskii (KK). 
Recently we have shown that the KK Hamiltonian does not support long range spin order at any 
nonzero temperature due to several novel hidden symmetries that it possesses. Here we extend our 
theory to show that these symmetries also apply to the underlying three-band Hubbard model. 
Using these symmetries we develop a rigorous Mermin- Wagner construction, which shows that the 
three-band Hubbard model does not support spontaneous long-range spin order at any nonzero 
temperature and at any order in t/U ~ despite the three-dimensional lattice structure. Introduction 
of spin-orbit coupling does allow spin ordering, but even then the excitation spectrum is gapless due 
to a subtle continuous symmetry. Finally we showed that these hidden symmetries dramatically 
simplify the numerical exact diagonalization studies of finite clusters. 



PACS numbers: 75.10.-b,71.27.+a,75.30.Et,75.30.Gw 



I. INTRODUCTION 



The transition metal oxides have been the source of many fascinating physical phenom- 
ena such as high temperature superconductivity^ , colossal magnetoresistance^, half metallic 
perovskites^ and orbiton physic3^>^. These surprising and diverse physical properties arise 
from strong correlation effects in the 3d bands. As a first-step towards a better understand- 
ing of these systems, there is much recent interest in the magnetic properties of transition 
metal oxides with cubic ABO3 structure, where three (i-orbitals are degenerate^ (see inset 
to Fig. 1). In cubic oxide perovskites, the crystal field of the surrounding oxygen octa- 
hedra splits the (i-orbitals into a two-fold degenerate Cg and a three-fold degenerate t2g 
manifold. In most cases, these degeneracies are further lifted by a cooperative Jahn- Teller 
(JT) distortion^, and the low energy physics is well described by an effective superexchange 
spin-only modei^^. However, some cubic perovskites, such as LaTiOs, do not undergo a 
significant JT distortion, in spite of the orbital degeneracy^. We will mainly consider the 
simplest "idealized" version of this model, in which the magnetic ions occupy sites whose 
local site symmetry is at least tetragonal. In these systems, the effective superexchange 
model must deal with not only the spin degrees of freedom but also the degenerate orbital 
degrees of freedo m^i^i^^ . The large degeneracy of the resulting ground states may then yield 
rich phase diagrams, with exotic types of order, involving a strong interplay between the 
spin and orbital sectors^»i2*ii. 

In the titanates, there is one d electron in the t2g degenerate manifold, which contains 
the wavefunctions |x) = dyz, \y) = d^^, and \z) = d^y. We will refer to |?/), and \z) as 
"flavors." Following Kugel and Khomskii (KK)Ji, one starts from a Hubbard model with on- 
site Coulomb energy U and nearest-neighbor (nn) hopping energy t. For large [/, this model 
can be reduced to an effective superexchange model, which involves only nn spin and orbital 
coupling, with energies of order e = /U . This low energy model has been the basis for 
theoretical studies of the titanates . ^i^?'^^'^^ In particular, it has been suggested^^ that the KK 
Hamiltonian gives rise to an ordered isotropic spin phase, and that an energy gap in the spin 
excitations can be caused by spin-orbit interactions^^. However, these papers are based on 
assumptions and approximations which are hard to assess. Recently^ (this will be referred 
to as I) we have presented rigorous arguments which show several unusual symmetries of 
the KK Hamiltonian. Perhaps the most striking symmetry is the rotational invariance of 
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the total spin of the electrons in a-flavor orbitals in each plane perpendicular to the a-axis. 
As reported previously in I, this symmetry implies that the system does not support long- 
range spin order at any nonzero temperature, despite the underlying three-dimensional lattice 
structure. Inclusion of spin-orbit interaction destroys the independent rotational invariance 
of the spin associated with each orbital. Although long-range order at nonzero temperature 
occurs when this perturbation is included, the spin system still possesses enough symmetry 
that the excitation spectrum remains gapless.i^*2iil (This conclusion might be surprising, 
because once spin-orbit interactions are included, the system is expected to distinguish 
directions relative to those defined by the lattice.) The purpose of the present paper is to 
extend the theory for the KK model presented in I to the underlying three-band t2g Hubbard 
model, and to elaborate further on the consequences of these results. 

Briefiy this paper is organized as follows. In Sec. II we introduce the Hubbard model, 
which forms the starting point for "idealized" treatments of these systems, and we discuss 
the KK Hamiltonian which is its byproduct. In Sec. HI we derive the symmetry results 
both on the original Hubbard model and also on the KK Hamiltonian, which are central to 
most of our arguments. In Sec. IV we give the details of the Mermin- Wagner construction, 
which enables us to rigorously conclude that this model does not support long-range order. 
In Sec. V we show how these unusual symmetries lead to a dramatic simplification in the 
numerical determination of the ground state of a cube of eight sites governed by the KK 
Hamiltonian. In Sec. VI we give a symmetry analysis which shows that even in the presence 
of spin-orbit interactions the excitation spectrum must have a gapless Goldstone mode. In 
Sec. VII we summarize our work. Our major conclusion is that to understand real physical 
systems, such as LaTiOa, which do show long-range spin order and which do have a gap in 
their elementary excitation spectrum, it is crucial to include perturbations which destroy 
the symmetries described here. 

II. THE HAMILTONIAN 

The system we treat is a simple cubic lattice of ions with one electron per ion in a d 
band. Following the seminal work of Kugel and Khomskii— (KK) , we describe this system 
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by a Hubbard Hamiltonian Tin of the form 

ia ij a/Bcr i a(3 

where c|„^ creates an electron in the orbital labeled a in spin state a on site i and Na{i) = 
J^a^laa'^iaa, = — {Uaaf^), whcrc is the Crystal field energy of the a orbital, tajsii^j) 
(which we assume to be real) is the matrix element for hopping between orbital a of site i and 
orbital (3 of site j, and Uap is the direct Coulomb interaction between electrons in orbitals a 
and (3 on the same site. Note that this Hamiltonian does not include the somewhat smaller 
Coulomb exchange terms . [We write Tin in the above form to emphasize that apart from 
the hopping term, the Hamiltonian is a function of Na{i).] We assume that the crystal field 
splits the five orbital d states into three low-energy states, dyz = d^z = \y)i and 
dxy = \z) and that the two other Cg states have high enough energy that we can neglect 
their presence (Fig. 1). This structure is consistent with tetragonal or higher site symmetry. 
(For strictly cubic site symmetry the higher two states form a degenerate doublet and the 
lower three states a degenerate triplet. For lower-than-cubic symmetry the degeneracy in 
energy of the states is removed but, as long as the symmetry is tetragonal or higher, the 
wave functions are those listed.) We confine our attention to the wide class of materials 
in which hopping between magnetic ions is mediated by intervening oxygen ions. In that 
case, following KK, we note that the hopping matrix taf3 between nn's is diagonal in orbital 
indices and also t^a = if the nn bond lies along the a axis, which has been callec^ the 
"inactive" axis for a hopping. Later on we will discuss modifications caused by further- 
than-nearest-neighbor hopping via oxygen ions. The symmetry of the nn hopping matrix 
is illustrated in Fig. |2l Note that these symmetries only hold for tetragonal or cubic site 
symmetry. They are broken by a rotation of the oxygen octahedra. 

When t <^U , KK reduced the above Hubbard Hamiltonian at lowest order in t/?7 to an 
effective Hamiltonian for the manifold of states which remain when U oo. We will call 
this Hamiltonian the KK Hamiltonian, and it will be denoted T^kk- This Hamiltonian can be 
regarded as a many-band generalization of the Heisenberg Hamiltonian, which is obtained 
from the single-band Hubbard model with one electron per site, for which the exchange 
constant is J = At'^/U. To make the analogy with the Heisenberg model more apparent, the 
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KK Hamiltonian is often written in terms of spin variables as 

(ij) 

where Sj = ^cTj is the vector spin operator for an electron on site i, {ij) indicates that the 
sum is over pairs of nearest neighbor ions, e = t'^/U, and the exchange "constant" is now an 
orbital operator written in terms of spinless fermion operators as 

Jij = X/ X/ (4a(^il3^][3^ja , (3) 

where al^ creates an electron in orbital |a) on site i (of either spin) and a 7^ ij means that 
we sum a over the values x, y, and z, except for the value of a which corresponds to the 
coordinate direction of the bond ij, which is the inactive axis for orbital a. 
For our purposes it is more convenient to write Hkk = 'Hx + 'Hy + 7i^, where 

where {ij) G a means that the sum over pairs of nearest neighbors is restricted to those for 
which Yij is along an a-axis. It should be clear that TYkk will inherit whatever symmetries 
which are present in T^h, although additional symmetries are to be expected. For instance, 
one can regard the elimination of states in which any site is doubly occupied as resulting 
from a canonical transformation which eliminates such states.^" Such a programme can be 
carried out for the single band Hubbard model, even when it is not at half filling; The 
important point is that whatever rotational symmetries we uncover in Tin should apply 
also to T^KK- On the other hand, the conservation law that each site has a single electron 
will only hold for T^kk because it is precisely the double occupancy sites which have been 
eliminated by the canonical transformation. 

III. HIDDEN ROTATIONAL SYMMETRY 

Previously in I^^ we pointed out several unusual symmetries of the effective Hamiltonian 
TYkk- Here we extend our analysis to the underlying Hubbard model. It is useful to rewrite 
Tin to display explicitly the form of the hopping matrix element: 

= E + ^E' E E cLc,a. + ^ E E Uo.pN^{i)Np{i) , (5) 

iaa ij a^ij o" i af) 
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where the prime on the sum over i and j hmits this sum to the case when these sites are 
nearest neighbors of one another. 

First we show that the total number of electrons in a orbitals in any chosen plane of sites 
perpendicular to the a-axis is conserved. For that purpose we note that the only operator 
which changes the occupancy of orbital states is the hopping term, T, where 

T = ^E' E E 4aaCjaa • (6) 

Electrons in an a-orbital may hop only to another a-orbital. Thus, the total number of a- 
electrons (by this we mean electrons in a orbitals) is a good quantum number Furthermore, 
hopping between a-orbitals can only take place within the same a-plane, i.e. within the same 
plane perpendicular to the a-axis, which is the inactive axis for a-flavor electrons. Thus, the 
total number of a-electrons in each a-plane is a good quantum number, so that the operator 
A^a(ao) commutes with the Hamiltonian, where 

^"("o) = EE^aaQaa- (7) 

Here the notation z G ao indicates that i is summed over all sites in the a-plane for which 
Via = ao, where is the a component of the position of site i. As mentioned, this property 
must be inherited by the KK Hamiltonian, and indeed one can show this explicitly. 

We now consider the much more general symmetry induced by an arbitrary unitary 
transformation among spin states of a given orbital flavor. We consider the effect of the 
transformation applied to operators on site i, 



E (8) 

V 

and 

Ciacr = U^^CjQ,^ = \jli^Ciar) , (9) 

V n 

where U is an arbitrary two dimensional unitary matrix and a is a fixed flavor the choice 
of which is arbitrary. Note that 



cr 

= E [u^u 

VP 

= J2^lav^iav = Nai-l) ■ (10) 
V 



crrip 



This shows that the first and third terms of Eq. (jS)) are invariant under this type of local 
unitary transformation applied only to orbital a on site i. The second term will likewise 
be invariant under an arbitrary U2 transformation, providing we transform all operators 
connected by hopping in the same way. This means that we set 

4/3. = Eu..(^,/3)£U (11) 

and 

Ci/3a = I] Uj„ Ci/3r, , (12) 

n 

where U(2,/5) is the unit (2 x 2) matrix, unless both /? = a and i is in the given a-plane: 
"^ia = ttO) in which case U(i,/3) = U, where U is an arbitrary two dimensional unitary 
matrix. Thus we apply the transformation U to all electrons of a given orbital flavor which 
are in an arbitrarily chosen plane perpendicular to the inactive axis for this flavor. To see 
what this means, we write the spin operator for a-electrons in the chosen plane as 

Sq!(«o) = c]ni^(T„j^CiQ,ri = ^ia ; (l^) 

where Sia is the vector spin operator for electrons on site i of orbital flavor a. Then we 
conclude that since U is an arbitrarily chosen unitary matrix the above discussion shows 
that Sa(Q;o), the total spin summed over all electrons in orbital a which are in any arbitrarily 
chosen a plane can be rotated arbitrarily at zero cost in energy. In other words the total 
spin of this orbital flavor in an arbitrarily chosen plane perpendicular to the inactive axis of 
this flavor, as well as the ^-component of this total spin, are also good quantum numbers. 
Applying the transformation PHI and (fT^ to the KK Hamiltonian (j3)) we note that, as 
should be the case, this Hamiltonian is also invariant under the rotation of the total spin of 
all electrons in orbital a in a given a-plane. 

We now discuss these results. First of all, the KK Hamiltonian of Eq. Q is somewhat 
similar to the Heisenberg Hamiltonian in that it involves the scalar product of spin operators 
on different sites. This form guarantees that the Hamiltonian is invariant under a global 
rotation of spin. We have demonstrated a much stronger symmetry, which only holds because 
the constant added to <Tj ■ (Tj in Eq. Q is unity. [Only in this case is it possible to write Eq. 
(0} in the alternate form of Eq. Q.] Furthermore, had we allowed Coulomb exchange in the 
Hubbard Hamiltonian, that is, had we allowed spin exchange in the Coulomb interaction 
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term of Eq. (0), then we would have global spin rotation (just as in the Heisenberg case), 
but we would not have had the invariance against independent rotations of the spins of each 
flavor. In fact, any interaction which allows an electron of one orbital flavor to convert (via 
hopping or some other interaction) into a different orbital flavor will also clearly invalidate 
the property that allows independent spin rotations for different orbital flavors. Probably 
the most important perturbation which allows such off-diagonal hopping in orbital flavor is 
the non-coUinear M — O — M bonding due to the rotation of the MOg-octahedra in real 
transition metal oxides. This rotation angle is 156° in LaTiO^^ and therefore it has to be 
taken into account in any theory to explain the observed experimental properties of LaTiOs. 
Similarly, this symmetry is also destroyed by spin-orbit interactions. Coulomb exchange 
interactions, hopping between second nearest neighbor oxygen ions, and less importantly by 
direct exchange between nearest neighboring Ti ions and by dipolar interactions. 

IV. MERMIN- WAGNER CONSTRUCTION 

In I we used a Mermin- Wagner construction^^ to establish the absence of long-range 
spin order due to the KK Hamiltonian in three (or less) dimensions, but we did not give 
full details of this construction. Here we present such a construction, but obtain a more 
powerful result by working with the Hubbard Hamiltonian. Results obtained from the KK 
Hamiltonian are only valid to leading order in t/U . By working with the Hubbard model 
Hamiltonian, we will establish these results to all orders in t/U for the Hamiltonian of Eq. 

The Mermin- Wagner construction relies on the Bogoliubov inequality,— which is 

^([[C,7^]_,Ct]_)({AAt}+) > A:T|([C,A]_)p . (14) 

where (X) denotes the canonical average of the operator X at temperature T, [X, y]_ is 
the commutator of X and Y and {X, y }+ is the anticommutator of X and Y . As we shall 
see, it is also crucial to note that 

Cn^{[[C,'HUC^]-)>Q , (15) 

unless [C, 7i]_ = 0, in which Ct-c = 0. To verify this we write out the expression for Cn in 
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terms of the system's exact energy eigenstates, |n): 



n 



(16) 



{n\C\m){E^ - En){m\C^\n) - {n\C^\m) {m\C\n){En - 

with pn = e~^^'' /Z, where Z is the partition function, 

Z = 5]e-^^", (17) 



n 



and (3 = l/{kT). In the second term of Eq. (jl6j) we interchange the roles of n and m, after 
which we obtain 

Cn = J2\i^\CH\^[Pn-Pm][Em-Er,] . (18) 

n,m 

We have that 

[Pn-Pm][Em~En] = Z-h-f''~^"+^-^/^2En,mSmh[{(3En,m/2]) , (19) 

where En,m = En — Em- Thus the summand in Eq. (fTHjl is nonnegative. In fact, for Cn to 
be zero, one has to have {n\C\m) = whenever En 7^ Em- In other words, C-h is real and 
nonnegative and can only be zero if the operator C commutes with the Hamiltonian, a case 
we will not encounter here. 

For the purposes of this construction we will add to the Hamiltonian a term conjugate to 
the magnetization at wavevector K of electrons in a-orbitals, so that we treat the Hamilto- 
nian 

^ R R 

= nu- hNS^„ , (20) 

where Sr^ is the spin operator for the a orbital of site R, as in Eq. (fT^ . We will set the 
order-parameter wavevector K to be zero to exclude ferromagnetic long-range order and to 
be (vr, TT, 7r)/a to exclude antiferromagnetic long-range order, where a is the lattice constant. 
However, as can be seen later, we may choose K arbitrarily, in which case our construction 
excludes the possibility of long-range order at whatever K is chosen. 
To rule out long-range spin ordering, we apply Eq. (fT^ taking 



a = Ee-^'^-^'cLrCRai , (21) 
R 



= Ee'^''^''^-^cl^iCna, , (22) 

R 



where R is a lattice site. Note that the operator Cj^q,|Cr,q,; can be identified as S^^, the 
raising operator for Sr^- Then Eq. (|T^ gives 

4}.) > 2kT E Jfi^'lf '^^^If , • (23) 

Here and below the sum over k is over the wavevectors of the first Brillouin zone, where 
is the total number of sites. We have 

k k R S 

The sum over k gives NSn^s, so that 

E({^k,4}+) = A^E({cLiCR„T,cLTCRai}+) • (25) 

k R 

Here the sum over R consists of a sum of 2N products, each of which is bounded by unity. 
Thus we have 

2iV^>E({^k,4}+) (26) 



which gives 



^2 > fcTV . (27) 



Now we evaluate 



R,S 

= E e^'' ''[cLTCR"i, cLiCR^t]- • (28) 

R 

Using 

[AB, CD]. = A{B, C}+D - {A, C}+BD + CA{B, D}+ - C{A, D}+B . (29) 
we then get 

[cRajCRa;, Cj^^^Cr„|]_ = Cj^at'^R"! " '^Rai'^R-"i (^O) 



10 



(which is nothing more than [S^, ^na]- — 2'S'rq,), so that 

([a, Ak]-) = 2^e^^-^.,(a,R) = 2NS^^ . (31) 

R 

For K — S'kq, is the ^-component of the spin per site of a-flavor orbitals and for K = 
(tt, tt, 7r)a it is the z-component of the staggered magnetization per site of a-flavor orbitals. 
(For general K it is the amplitude of the Fourier component of the z-component of spin of a- 
fiavor orbitals at wavevector K.) We concentrate on the ferromagnetic and antiferromagnetic 
cases. We have the inequahty 

l>AkT\SU'^{X^)-\ (32) 

k 

where is positive real and is given by 

= [[c^,nuclu . (33) 

Now we need the double commutator of Ck with 7i. The Hamiltonian is the sum of four 
terms: Ti involving the hopping matrix element tfj, T2 involving the crystal field energies e^, 
T3 the on-site Coulomb interaction, and T4 the interaction with the staggered field h > 0. 
Since T2 and T3 depend only on Na{i), one can establish that 

[Ck,r2]- = [Ck,T3]_ = 0. (34) 

Also 

[Ck,T4] 

so that 

{[[C^,T,]_,Cl].) = (/iEe^''-''[4aTCRai,cLiCaaT]-) 

R 

R 

= 2NhS^, . (36) 
The crucial calculation is the commutator involving Ti. 

Fk = [Ck, ri]_ = Yl e-'^'-^icL^cn^l: E T)cL.CTa.]- , (37) 

R S,T,CT 



'-RaT'^R-aJ- 



(35) 



R 
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where R, S, and T are lattice sites and we dropped terms involving orbitals other than a 
since they obviously commute with Ck- Also in this section we set ta(R, S) = ta(Il — S) to 
denote the hopping matrix element (assumed to be nonnegative for convenience) between a 
orbitals on sites R and S and we use tQ,(R) = tQ,(— R). 
We have 



RSTo 



RST 



1 -e 



where <5 is a nearest-neighbor vector. Then 



nS 



(38) 



= [[Ck, Ti]_, + 2NhS^^ = [Fk, 4]- + '^NhS- 



E' 

RS(5 



^ik-(S-R) 



1-e 



Ka 



E^4^) 

R^ 

E^"(^) 

R^ 



1 - e 



ik (5 



+ 2NhS^^ . (39) 



According to Eq. (fTH|) the quantity (Xk) is real and positive (no matter what the value 
of k) and therefore we can maintain the inequality (j32p if we replace the sum S, where 

1 

k W ' 



(40) 



by a quantity which is less than S. To do that we write 

S = ^ I ^ 

2 ^Xk) ^ (X_k) 

9 

> E 



(41) 



^ (Xk) + (x_k) • 

This follows from the fact that when a and b are positive, then + b^^ > 4/(a + b) 
Thus 



l>8A:T|5^j2E[W + (X_kr' 



(42) 
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Here 



[(Xk) + (x_k)] /2 = j:us) 



R<5 

+2NhS^^ 



[cos(k ■ (5) - 1] ((cLtCr+5„^) + (cLiCR+<5ai) 



(43) 



Since the quantity on the right-hand side of Eq. ()43p is positive, we may replace it 
(without increasing S) by the sum of the absolute values of bounds to its terms, so that 
finally Eq. ^ yields 



-1. 



^ \ s 



t/I 



(44) 



where 



N 



He^MIi- cosik-s)] + h\su] 

k \ A / 



-1 



(45) 



and we have used the bound 



where Co = 1. In I we gave the result obtained in an entirely analogous fashion from the 
KK Hamiltonian. That result is valid for the KK Hamiltonian, which itself assumes the 
validity of the expansion in powers of t/U. That result could be obtained here by noting 
that to lowest order in perturbation theory in t/U, the left-hand side of Eq. is given 
by24 Cq ~ 2\t\/U. The present result avoids any rigorous discussion of the validity of the 
expansion in powers of t/U, but uses the much cruder bound of Eq. pHjl with Cq = 1. In 
any event, an entirely analogous construction can be carried out for the KK Hamiltonian, 
as discussed in I. 

The analysis can now be less formal. We estimate the quantity / on the right-hand side 
of Eq. ()45|1 . For a = z, for example, we have 

^ta((5)[l - cos(k ■ d] = 2t[2 - cos(/c^a) - cos{kya)] . (47) 
S 

Since the sum over k is dominated by small |k|, we may write 

, 2k\dk\ 



k\^Bh\SUIt 



A In 



1 + 



Bh\S^ 



(48) 



K2 
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where A,B, and xq are constants of order unity, whose exact values need not concern us. In 
all we obtain the bound 

\S^,\ < A'T-^/^\\n\h\/t\-'/\ (49) 

where A' is a constant. Obviously, this bound applies also to any of the other components, 
S'kc, of the spin of any flavor a. This bound implies that as h ^ 0, we must have \S^^\ 0, 
for all choices of K, a, and fi. Thus we conclude that this Hubbard model (and perforce also 
the KK Hamiltonian obtained from it) can not support spontaneous long-range spin order at 
any nonzero temperature. Our rigorous arguments have nothing to say about orbital order. 
However, we believe that there is no spontaneous long-range orbital order for the KK model 
with no spin-orbit interactions. 

From the form of the bound we also expect that the lower critical dimension for the ap- 
pearance of long-range spin order is (i< = 3. As will be discussed elsewhere,— this conclusion 
can also be understood within a renormalization group analysis. 

The above development allows us to make some comments on whether or not there can 
be spontaneous breaking of parity symmetry. Note that the average of the quantity on the 
right-hand side of Eq. ()39|) must be positive. However, apart from the term proportional to 
h, this quantity vanishes for k = and has a term linear in k. The term linear in k gives a 
contribution to {X^) of 

= ^E^"(^)(k-^)(cLTCR+5„t -cLiCR+<5^^) . (50) 

R(5 

This quantity must be real, but must also go to zero as /i — > 0, in order that (X^) always 
be positive. If parity is preserved, then, of course, {cn.aa'^ii+Saa'^ ~ ('^Rao-'^R-^ao-) 
term linear in k does vanish. The converse, is not quite proven, because the quantity in 
Eq. ()50|) could vanish without parity being maintained. So, this development rules out a 
spin-independent breaking of parity. In addition, we remind the reader that we can choose 
K arbitrarily, after which the above argument rules out spontaneous helical spin order. 

V. EXACT DIAGONALIZATION 

In this section we show that the hidden symmetries discussed above are very useful 
in simplifying the exact numerical studies of small clusters. Such cluster studies are very 
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important, because they provide a way of checking analytical results and also give insight into 
the nature of the ground state when there is no long range order. The main problem in such 
studies of finite clusters is usually the large matrix sizes that have to be diagonalized. For 
example, to treat the KK Hamiltonian for a cube of eight sites requires the diagonalization 
of a matrix of dimensionality 6® ~ 1.7 million. Since the KK Hamiltonian is rotationally 
invariant, one can obtain the full spectrum by working within the subspace of Sf = 0, 
where the size of the Hamiltonian matrix is now reduced to 70 x 3^ ~ 1/2 million. The spin 
degeneracy of each individual state in this manifold is 2S + 1, where S{S + 1) is the square 
of the total spin of the wavefunction. 

Initially we obtained the low-energy spectrum of eigenstates for a cube of eight Ti ions 
by diagonalizing the 1/2 million dimensional matrix described above. Since the Hamilto- 
nian matrix is very sparse, it is possible to obtain the eigenvalues and the corresponding 
eigenvectors in a small interval of the spectrum starting from the lowest eigenvalue, using 
a standard sparse matrix diagonalization routine. The ground state is found to be three- 
fold degenerate, with energy — 6.2716e. By analyzing the wavefunctions we found that the 
ground state had a total spin S = and either A^^^ = 4, A^y = 4, A^^ = 0, or the two cyclic 
permutations of these quantum numbers, where A^^ is the total number of a-flavor electrons. 

The same energy spectrum and the corresponding wavefunctions can also be found by 
diagonalizing much smaller matrices if one works within a manifold defined by the conserved 
numbers applied to each face of the cube (see Fig. E)). We actually have eighteen conserved 
quantum numbers and there is an astonishing numerical simplification when maximal use is 
made of these symmetries. For example, the ground state found by dealing with the matrix 
of dimensionality 1/2 million could alternatively be found by diagonalizing the Hamiltonian 
matrix within a manifold of just 16 states! Given the fact that the ground state has S = and 

= Ny = 4, one can easily conclude that this corresponds to having A^^, i = 2 electrons in x 
orbitals in the first x-plane and Nx^2 = 2 electrons in x orbitals in the second x-plane. For any 
asymmetric choice, such as A^^^^i = 1 and Nx^2 = 3, for instance, one would have an additional 
degeneracy associated with interchanging N^^i and Nx,2- Thus we were sure that these states 
found numerically had to have the quantum numbers N^^i = A^a;,2 = = ^y,2 = 2. In 
addition, the total spin of the two x electrons in each x-plane had to be zero in order to be 
consistent with the lack of spin degeneracy. Thus it was clear that the ground state had to 
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consist of a sum of terms, each term being a product of four spin singlets. We define 

^ 2-i/^[cLt4.^ - cLi4^]|vac) = , (51) 

with a similar definition of {ij)y, where |vac) is the vacuum state. Then the ground state 
consists of a linear combination of terms, each term being a product of four dimer singlets 
which we write as 

\{ij)x{kl)^{mn)y{op)y) , (52) 

where the sites i and j are in the first x-plane, k and I are in the second x-plane, m and n are 
in the first y-plane, and o and p are in the second y-plane. The dominant configurations in 
the ground state are the dimer states with the lowest expectation value of the Hamiltonian, 
namely the states in panels (b) and (c) of Fig. El From one of these, say, the one shown in 
panel (c), one can generate the manifold of states which can be obtained by hopping along 
one or more of the z-directed bonds, labeled a, b, c, and d in the figure. The Hamiltonian 
matrix within this manifold of 16 states, is given explicitly in Table I.i^ To construct this 
matrix we used the fact that 

malKKKma) (53) 

is — e if sites i and j are nn's along an axis which is not inactive for a orbitals. Otherwise 
this matrix element is zero. (For instance, diagonal dimers like the X and Y dimers involved 
in the hopping shown in panel (d) of Fig. 01 make zero contribution to the diagonal matrix 
element. So configuration (d) has a diagonal matrix element of — 2e due to the nn X and Y 
dimers.) Also the matrix element 

{{ij)ximn)y\n\{mj)^^{in)y) (54) 

has the value — e if the sites i and m which are involved in the hopping are nn's along the 
active (z) direction, and is zero otherwise. It is amazing to us that the original Hamilto- 
nian, which was obtained from a matrix of dimensionality approximately 1/2 million, could 
be reduced by symmetry considerations to an eigenvalue problem of dimension 16. The 
existence of these symmetries was numerically confirmed when we found the same value 
for the ground state energy from the 16 dimensional matrix as from the full matrix. The 
ground-state wavefunction, in terms of the states in the order used for the matrix, is 

(a, /?, /?, /3, /3, 7, 7, (5, 5, 7, 7, /?, /3, /?, /?, a) , (55) 
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where a = 0.40915, (3 = 0.23235, 7 = 0.21758, and S = 0.14819. 



VI. SYMMETRY, SPIN-ORBIT INTERACTION, AND THE EXCITATION 
SPECTRUM 

Here we consider the addition of a perturbation which breaks the symmetry whereby each 
orbital flavor of spin can independently be rotated at zero cost in energy. For concreteness 
we consider the spin-orbit interaction. Possibly one's first intuition about the effect of spin- 
orbit interaction would be that it would cause the system to exhibit long-range order (this is 
correct) and that the elementary excitation spectrum would have a gap, because the orbits 
would define a set of favored directions. In a study of a similar Hubbard model for the 
cuprates this second conclusion was shown to be false.— In fact, there it was shown that 
addition of spin-orbit interactions to an isotropic Hubbard model did not lead to a gap in 
the spin- wave spectrum, but that a gap does result with the addition of both spin-orbit and 
Coulomb exchange interactions (these are sometimes called Hund's rule coupling). Here we 
will establish the first part of this scenario, namely that adding only spin-orbit interactions 
to the KK Hamiltonian or to the Hubbard model of Eq. (jT)) does not lead to a gap in the 
excitation spectrum. 

Following Ref. (Si we now introduce a transformation to pseudo-spin which yields a rota- 
tionally invariant Hamiltonian. We write^i 

V 

where dj^^ creates an electron in orbital a of site i with pseudospin i]. Since pseudospin is not 
the most intuitive concept, we will here give a discussion that avoids use of this quantity. In 
terms of "real" spin we will show that the Hamiltonian is invariant under the transformation 

T 

where 

V(") = a„Ua„ , (58) 
with an arbitrary unitary matrix, U. 
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We now consider the effect of this transformation on the Hubbard model of Eq. (P). 
The Coulomb interaction and single-site crystal-field energy are both clearly invariant un- 
der this transformation, because they only depend on the total number of electrons in each 
orbital, J2a clacr<^ioia, aud this quantity is not changed by this transformation. Since hopping 
is diagonal in the orbital indices, the hopping term is also invariant under this transforma- 
tion. Finally, we study how the additional spin-orbit interaction transforms. We write this 
interaction as 



^s-o = AEEE(«l^7l/5)[ 



(59) 



where {a\L^\i3) is the orbital angular momentum matrix element, and is nonzero only when 
a, (3 and 7 are all different: 



(a|L^|7) = -ie, 



a/37 



(60) 



Thus 



where 



= E("l^7l/^)MprcLpCi/3^ , 



M 



pr 



E[ 



HP 



pr 



pr 



(61) 



(62) 



The orbital matrix element (a|L^|/?) guarantees that all three of these indices are different. 

In that case we set 



(63) 



where X is the unit matrix and eap-^ is the antisymmetric tensor. Thus 



M 



pr 



pr 



(64) 
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Since this matrix has the same form as before the transformation, we have shown that the 
Hubbard Hamiltonian with spin-orbit interactions is invariant against the transformation by 
the unitary matrix of Eq. (j58|) . where U is an arbitrary unitary matrix. Thus the invariance 
with respect to arbitrary rotation of pseudospin gives rise to an arbitrary transformation of 
real spin. This indicates that there is a continuously degenerate manifold for the ground 
state, and therefore we do not expect a gap in the excitation spectrum. A rigorous analysis 
of the structure of the Goldstone modes is more complicated^^ and may require a statement 
of whether or not there is long-range order in the orbital sector, in order to say whether these 
modes are propagating or diffusive. Intuitively it seems likely that if there is long-range spin 
order, there should be propagating spin-wave modes, at least one of which by our argument 
will not have an energy gap at zero wavevector. 

One final point is worth noting. We have discussed that we expect long-range antifer- 
romagnetic spin ordering in the presence of spin-orbit coupling. Once we have antiferro- 
magnetic spin order, the spin-orbit perturbation acts like a staggered field on the orbital 
variable, thus inducing long-range antiferromagnetic order in the orbital variable (L). This 
type of orbital ordering is not the same as that in which the thermal expectation value of 
the orbital occupation numbers, (A^^^), become unequal to one another. 

VII. CONCLUSIONS 

In this paper we have analyzed several unusual symmetries of the three-band t2g Hubbard 
model and the associated low-energy KK Hamiltonian for three-fold degenerate orbitals in 
the ideal ABO3 cubic structure. As pointed out previously^, it is evident that the total num- 
ber of electrons in any of the three orbital states is a good quantum number. In addition,— 
the number of electrons in any arbitrarily chosen plane perpendicular to the a-axis which are 
in a orbitals is also a good quantum number. We showed that the Hamiltonian is invariant 
with respect to independently rotating the spins of any single orbital state a for all spins 
in an arbitrarily chosen plane perpendicular to the a-axis (which is the inactive axis for 
such orbital spins). These symmetries lead to a dramatic simplification in numerical diago- 
nalizations, as we illustrated by discussing exact diagonalizations for a cube of eight sites. 
Furthermore, using these symmetries we develop a rigorous Mermin- Wagner construction, 
which shows that the Hubbard model (and perforce the KK Hamiltonian derived from it) 
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does not support spontaneous long-range spin order at any nonzero temperature. These 
unusual symmetries are destroyed by almost any realistic perturbation. In the presence of 
spin-orbit coupling, for instance, we show that a special invariance with respect to rotation 
of pseudo-spin remains unbroken and due to this continuous symmetry we expect the excita- 
tion spectrum to have at least one gapless Goldstone mode. In analogy with the situation in 
the cuprates^ii^, the gap one might have expected in the presence of spin-orbit interactions 
will only develop in the presence of both spin-orbit and Coulomb exchange interactions.— 

Of course, in real systems, such as LaTiOs, long-range spin order and a gap in the 
elementary excitation spectrum are observed.— However, any credible expression for Tc, for 
instance, must involve a perturbation beyond the isotropic "bare" KK model, because, as we 
have shown, this model does not support long-range spin order. Elsewhere^^ we discuss the 
type of long-range order which mean-field theory predicts when various perturbations such as 
a) spin-orbit interactions, b) nnn hopping, and c) Hund's rule coupling are added to the KK 
Hamiltonian. In addition, an important factor to consider is the distortion and/or rotation 
of the oxygen octahedra surrounding each Ti ion. Similarly, to be credible any theoretical 
expression for the gap in the excitation spectrum must involve perturbations beyond spin- 
orbit interactions, because we have shown that with only spin-orbit interactions added to the 
"bare" KK Hamiltonian, the excitation spectrum is gapless. Finally, we note that it would 
be very interesting to study an experimental system with as small as possible deviations 
from the "bare" KK model. Such a system would show quite exotic properties (such as 
peculiar two-dimensional correlations of spins associated with different orbital flavors). 
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TABLE I: Hamiltonian matrix for the ground manifold in units of e. The states are specified by 
hops relative to the state |0) shown in Fig. (So a indicates hopping along the z-directed bond 
labeled a, and similarly for b, c, and d.) We introduce the notation a = abcda. 








a 


b 


c 


d 


ab 


be 


bd 


ac 


ad 


cd 


d 


c 


b 


a 








-4 


-1 


-1 


-1 


-1 



































a 


-1 


-2 











-1 








-1 


-1 




















b 


-1 





-2 








-1 


-1 


-1 


























c 


-1 








-2 








-1 





-1 





-1 

















d 


-1 











-2 








-1 





-1 


-1 

















ab 





-1 


-1 








-2 

















-1 


-1 











be 








-1 


-1 








-2 














-1 








-1 





bd 








-1 





-1 























-1 





-1 





ac 





-1 





-1 























-1 





-1 








ad 





-1 








-1 














-2 








-1 


-1 








cd 











-1 


-1 

















-2 








-1 


-1 





d 

















-1 


-1 





-1 








-2 











-1 


c 

















-1 





-1 





-1 








-2 








-1 


b 


























-1 


-1 


-1 








-2 





-1 


a 




















-1 


-1 








-1 











-2 


-1 






































-1 


-1 


-1 


-1 


-4 



23 




24 




y 



FIG. 2: Symmetry of the hopping matrix element for cubic site symmetry, (a) The hopping matrix 
between different flavors is zero, (b) The z-axis is inactive for z orbitals. 
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FIG. 3: Spin and orbital configurations for a cube of eight sites. The thick hues indicate singlet 
spin states (dimers) and X and Y indicate the orbital states of the electrons, (a) Schematic 
illustration of the quantum numbers associated with occupation of a orbitals in an a-plane. (Similar 
spin quantum numbers are not shown.) (b) and (c) Dominant configurations in the ground state 
wavefunction. Hopping between x and y orbital states is only allowed along ^-direction bonds and 
these bonds are labeled in panel (c). (d) A subdominant configuration in the ground state which 
is obtained from (c) by allowing the interchange of two {X and Y) electrons along the z-axis bond 
labeled "d" while each electron retains its membership in its original spin singlets (even though 
the position of the electron has changed). If |0) denotes configuration (c), then configuration (b) 
is |0) = ahcd\{)) and (d) is d|0) in the notation of Table I. We also show (e) and (f) states with the 
same quantum numbers but which are not coupled to the manifold of 16 states we analyze. These 
two states are actually eigenstates of the Hamiltonian with energy — 4e. 
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